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We study the applicability of the mass insertion approximation (MIA) for cal- 
culations of neutral meson mixing when squark masses are not degenerate and, in 
particular, in models of alignment. We show that the MIA can give results that 
are much better than an order of magnitude estimate as long as the masses are 
not strongly hierarchical. We argue that, in an effective two-squark framework, 
fh q = {rh\ + fh\)/2 is the best choice for the MIA expansion point, rather than, for 
example, fhg = (fhf + fh\)/2. 

I. INTRODUCTION 

The mass insertion approximation (MIA) is often used to simplify expressions involv- 
ing supersymmetric contributions to flavour changing neutral current processes from loop 
diagrams [p], 0, [|. The simplification is achieved by the replacement of a sum over all pos- 
sible internal propagators and the appropriate mixing at the vertices, with a single (small) 
off-diagonal mass insertion in a basis where all gauge couplings are diagonal. The result- 
ing expressions are formulated in terms of parameters which can be estimated in various 
supersymmetric models. 

It may seem, naively, that the smallness of the off-diagonal mass-squared matrix element 
would justify the approximation. The true picture, however, is that these off-diagonal ele- 
ments are the product of mixing angles at the vertices and mass-squared differences between 
intermediate squarks. The MIA, on the other hand, is a Taylor expansion only with respect 
to the latter, namely, the mass-squared difference (we give an exact formulation of these 
statements in section |T|). A small off-diagonal element does not necessarily imply a small 
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mass difference. Instead, it may be related to small mixing angles. But then the validity of 
the MIA is questionable. 

This is exactly the situation in the framework of quark-squark alignment (QSA) models 0, 
In this class of supersymmetric models the squark masses-squared are all of the same 
order of magnitude, a free parameter denoted by m 2 , and the mass squared differences 
between them are also of the same order of magnitude, that is: |m 2 — m 2 |/m 2 = 0(1). 
Apriori, this is a problematic situation for using the MIA. Yet, it is frequently used in the 
literature. 1 

We therefore study the validity of the MIA in the context of QSA models. We confirm that 
the approximation is applicable and useful for such models. We also clarify the connection 
between the (unknown) details of the squark mass spectrum and the MIA parameters. 

The organization of this work is as follows: We formulate the details of the MIA in 
section |I[ The analysis and our results for non-degenerate squarks masses are presented in 
section |T|. 



II. FORMULATION OF THE MASS INSERTION APPROXIMATION. 

Let us first formulate the details of the MIA. To illustrate it in the context of QSA, we 
study a specific example: The supersymmetric contribution to neutral K meson mixing from 
gluino box diagrams with two intermediate squarks. The relevant diagrams are shown in 
figure 0. We focus on the following term arising from these diagrams ||: 

MS D C (ziZ'lX^),) Hml ml m)) . (1) 

Here C is a numerical factor, given in terms of the K meson parameters: 

C EE (4?T)2 a ' m ^KV ^ 5 , 4 x 1Q 4 MeV3) (2) 

^2 z 
Zf- are the quark-squark mixing angles, and J 4 is given by 

11— 9 

T/~2~2~2\_- L r/~2~2~2\ l z, ~2w~2~2~2\ to\ 



1 In fact, since QSA models allow estimates of \mf — mj\/m 2 but not of the individual masses, the MIA 
provides the best way to derive meaningful results. 
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FIG. 1: Gluino box diagrams contributing to K° - K° mixing, 
with fhg the gluino mass, to*, m,- the down squark masses and 
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The MIA is nothing more than a Taylor expansion of J4. We choose to expand the squark 
masses around some point rh q . Owing to the specific form of J4 we find 



J 4 (m 2 g ,m 2 ,m 2 ) = 
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where x = m 2 /fh 2 , and Am| = mf—rh 2 . The coefficient C m+n (x)/mg < ' 1+m+r ^ is the (m+n)'th 
derivative of J4 evaluated at m q , times the symmetry factor l/(m\n\). (The exact form of 
the coefficient is not important for our purpose.) 
Substituting (j^) in (|l|) we get 

m,n=0 L m 9 J 

Note that a sum over i and j is implied. If the MIA is valid, this expansion converges fast 
and we can keep only the lowest order terms. However, due to the unitarity of Z d , the terms 
with either m = or n = vanish. The first non-vanishing contribution, therefore, will be 
from the term with m — n — 1. This term, however, is special since we can write (again, a 
sum over j is implied): 

Z% (Am*) Zft = Z d 2j (ml + Am?) Z* = {Mj) 2l , (8) 

where Mj is the squark mass-squared matrix in the basis where quarks masses and gluino 
couplings are diagonal. Thus, if the MIA holds we can replace (|7]) with: 

We stress that a small (M|) 2 i is not enough, by itself, to justify the use of the MIA. The 
question of validity should be considered in the context of eq. (H). We note, however, that 
the quality of the MIA is not completely equivalent to the quality of the approximation that 
is obtained by keeping only the lowest terms in (|6]). The reason is that the zeroth order term 
(m = n = 0), the first order terms (m = 1, n = and m = 0, n = 1) and some of the second 
order terms (m = 2, n = and m — 0, n — 2), while appearing in (^), do not contribute to 
the mixing amplitude in (0) due to the unitarity of Z d . In other words, the validity of (||) 
has to do with one of the second order terms in (|6|), rather than with all lowest order terms. 

Nonetheless, it is obvious that when the expansion parameter is small, \/\m 2 /rh 2 \ 1, 
the approximation is good. It is the condition in QSA models, \Arh 2 /rh 2 \ ~ 1 which needs 
a special consideration. 

III. THE CASE OF NON-DEGENERATE MASSES 

In order to study quantitatively the non-degenerate case, we simplify the form of Z d 
by assuming that only one mixing angle, namely Zf 2 , is large. This is usually the case 
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in QSA models. Such an assumption allows us to consider only the first two generations. 
Expression (|1]) then simplifies to 

C X cos* Bin* X [M^, ih\,m\) + Urn], ffif, 5$ - J 4 (ffi», ffij, ffi») - J 4 (S*, m*, fBJ)] , 

(10) 

where sin# « Zf 2 is the single large mixing angle. Equivalently, (]?]) can be written as 

tfxcos'Wflx f^M_) x ((Amr-N) m )(N)"-N)1 • (11) 
\ m ? / 

The above expression manifestly demonstrates the vanishing of terms with either m = or 
n = 0. The MIA is obtained by keeping in (|TT|) only the term with m = n = 1 



C 2 (m 2 /m 2 )\ f m 
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We can now test the accuracy of the MIA by defining the deviation parameter (using the 
symmetry of J4 with respect to mf and in 2 ) 

C 2 (ih 2 Jih 2 q )\ fm 2 -m 2 ^ 2 



m 2 J \ m 2 



r = 1 J: S— v _ 1 / —— _ _ . (13) 

Ji{m 2 , m\ , m\ ) + J^{m 2 , m 2 , m|) — 2 J^{m 2 , m\ , m 2 ) 

Figure shows function of mi and m 2 for m q = 2 TeV, m g = 1 TeV. 

Obviously, when rh\ pa ih 2 pa ih q , we expect the deviation to be small. The interesting 
result, however, is that the deviation is small (of order 10%) even for non-degenerate rh\ and 
ih 2 , as long as (mi + m 2 )/2 pa m g . This is the case, for example, when fh\ = 1 TeV, ih 2 = 3 
TeV and ih q = 2 TeV. In other words, the result of the exact expression (|10|) using ihi and 
m 2 , can be reproduced using the MIA of (|T2| ) with ih q pa (mj + m 2 )/2 and the knowledge 



of M 12 ~ cos 6 1 sin 6(ihf — m^/ih 2 . 

We see, therefore, that as long as the squark masses are not strongly hierarchical, the 
'right' choice of the expansion point ih q results with the MIA being a rather accurate approx- 
imation. This is very useful in QSA models where we have (m 2 —ih 2 ) /in 2 = 0(1). The ampli- 
tudes in these models can therefore be expressed in the MIA by using ih q pa (mi+m 2 )/2 ~ in. 
Although this result was demonstrated here for a specific mixing contribution and using spe- 
cific assumptions on the mixing angles, we found it to be quite a general result. 2 



2 Moreover, the example given here directly applies to the most stringent test of QSA models, namely, the 
contribution to D° - D° mixing. 
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FIG. 2: r as a function of m\ and m,2 for m q = 2 TeV, m g = 1 TeV. 

The result rh q ~ (mi + m-i)j c l is not trivial. Looking at (|TT|), it seems that the best 
strategy is to choose in 2 , = (rh\ + m 2 )/2. This choice results with Arhf = —Am 2 , which 
eliminates all terms with either m or n even. 3 In particular, it eliminates the next-to- leading 
order terms with either m = 1, n = 2 or m = 2, n = 1. Naively, one would expect faster 
convergence in this case and therefore obtaining a good approximation. 

This naive expectation is, however, wrong. Studying the form of the coefficients C m+n 
we find that the sign flips between even and odd m + n terms. The choice m 2 = (rh\ + 
m^)/2, which eliminates all odd m + n terms, eliminates therefore all the negative sign terms 
in (|TT1) . Since the C m+n coefficients decrease slowly, this induces a larger error and a worse 
approximation. 

On the other hand, we find that the choice m q ~ (mi + rri2)/2 is the most sensible one 
since it leads to an approximate cancellation between the next-to-leading order and the 
next-to-next-to-leading order terms. The extent to which such a choice is optimal, as can 
be seen in figure 0, is remarkable. 

Although we presented here explicitly only the LL and RR contributions to the mixing 

3 In other words, all terms with m + n odd are eliminated, since m + n odd implies either m even and n 
odd or vice versa. Some of the m + n even terms which are due to both m and n even are eliminated as 
well. 



7 



amplitude, the approximation also holds (over a somewhat smaller range of masses) for the 
LR and RL contributions. 

To summarize, we confirm that the common practice of using the MIA in supersymmetric 
models is justified even in the case of non-degenerate squark masses, as long as the masses 
are not strongly hierarchical. We showed that over a wide range of masses, the best strategy 
is to choose the MIA expansion point rh q to be the average of the masses involved. Using this 
strategy, the MIA provides a surprisingly good calculation of neutral meson mixing in models 
where a single mixing angle dominates, such as QSA 0, f| and 'effective supersymmetry' |7], 
1,1,0,0. 
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